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Abstract 

We consider the production, at future lepton colliders, of charged Higgs 
pairs in supersymmetric models. Assuming a relatively light SUSY sce- 
nario, and working in the MSSM, we show that, for cm. energies in the 
one TeV range, a one-loop logarithmic Sudakov expansion that includes an 
"effective" next-to subleading order term is adequate to the expected level 
of experimental accuracy. We consider then the coefficient of the linear 
(subleading) SUSY Sudakov logarithm and the SUSY next to subleading 
term of the expansion and show that their dependence on the supersym- 
metric parameters of the model is drastically different. In particular the 
coefficient of the SUSY logarithm is only dependent on tan (3 while the next 
to subleading term depends on a larger set of SUSY parameters. This would 
allow to extract from the data separate informations and tests of the model. 
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I. INTRODUCTION 



In the last few years, a considerable effort has been devoted to the precise formulation 
of the theoretical predictions for electroweak effects in pair production at future lepton 
colliders. In particular, the considered center of mass (cm.) energies have been those 
that represent the final goal of two proposed future machines, roughly one TeV for LC 
[I] and three TeV for CLIC [2]. The main motivation of the various investigations has 
been the fact that, within the electroweak sector of the SM, for cm. energies of the 
few TeV size, it has been realized [3-5] that unexpectedly large virtual effects arise at 
the one loop level, that could make the validity of this (relatively) simple perturbative 
calculation highly debatable. These terms have the analogous dependence on energy as 
those originally determined in QED by Sudakov [6]; at one loop, they can either be of 
squared logarithmic (leading) (DL) or of linear logarithmic (subleading) (SL) kind, and 
their numerical effect in several observables breaks the "safety" few percent limit (fixed 
by the aimed ( one percent) experimental accuracy) when one enters the few (2-3) TeV 
region, making the request of a higher order calculation to become imperative in that 
range. 

Within the SM framework a resummation to all orders actually exists to 
next-to subleading order for final massless fermion pairs, and to subleading order for gen- 
eral massive final pairs [7]. In the (particularly relevant) case of massive final fermions a 
comparison (to subleading logarithmic accuracy) between one-loop and resummed expan- 
sions has been performed [8]. The results indicate that, for cm. energies entering the few 
(2-3) TeV range, the discrepancies between the two approximations become intolerably 
(i.e. beyond a relative ten percent) large. On the contrary, for energies in the one TeV 
region, no appreciable difference shows up: to subleading order, the one-loop description 
appears there adequate. The question remains that of whether possible next-to sublead- 
ing (for instance constant) terms might play a role. 

For massive (third generation) final quark pairs, this problem was investigated in an 
"effective" way [9], trying to fit the exact one- loop calculation of a simple class of dia- 
grams with a logarithmic expansion also containing an additional constant. The result 
was that, in the one TeV region, this fit was adequately describing the exact calculation, 
with the logarithmic coefficients exactly predicted by the Sudakov expansion and a con- 
stant term given by the fit. The size of this term was (relatively) "large", about one half 
of the logarithmic contribution, with opposite sign (thus decreasing the overall effect). 
The apparent conclusion was that, in that energy range, a one- loop Sudakov expansion 
implemented by an extra constant term seems to be able to reproduce the exact calcu- 
lation with an accuracy that is largely sufficient, at the expected experimental precision 
level of a relative one percent. 

As a comment on the previous conclusion, it can be noticed that the possible exper- 
imental determination of the separate logarithmic and constant terms would not lead to 
any new information in the SM, since all the various coefficients will depend on the (fixed) 



2 



known values of the SM parameters. In this sense, accurate measurements of massive (and 
massless) fermion pair production at LC in the one TeV region can only provide, in the 
SM framework, another experimental test of the validity of the model. 

A natural question that arises at this stage is that of whether similar, or different, con- 
clusions can be drawn for the simplest supersymmetric generalization of the SM that can 
still be treated perturbatively, i.e. the Minimal Supersymmetric Standard Model. This 
model and its Sudakov expansion has been actually already considered to subleading log- 
arithmic accuracy in a number of papers, to the one-loop level for massless and massive 
fermion pairs [4,5] and resumming to all orders for sfermion and Higgs production [10]. 
In the latter case, a comparison between the one-loop and the resummed expansion has 
also been performed under the assumption of a relatively light SUSY scenario, showing 
that, in strict analogy with the SM situation, the two calculations are essentially identical 
in the one TeV region, while deep differences show up in the (2,3) TeV energy range. No 
effort was made to try to estimate the size of a possible next to subleading term in the 
one-loop expansion in the one TeV region from a fit to the exact calculation, analogous 
to that performed in [9] for the SM case. 

The aim of this paper is precisely that of studying the feasibility and the possible 
advantages of performing an effective logarithmic one-loop expansion, implemented by a 
next-to subleading term, in the energy region around 1 TeV, for the MSSM. We anticipate 
that a preliminary necessary condition will be that of a light SUSY scenario, in which 
all the SUSY masses relevant for the considered process are supposed to be not heavier 
than a few hundred GeV. This means that the starting picture will be one where SUSY 
has already been discovered via direct production of (at least some) sparticles, so that a 
number of SUSY masses is already known and fixed. This is not obviously true for other 
quantities of the model, like for instance tan (3 or other parameters of the Higgs sector. 

The purpose of our investigation will be that of showing that one fundamental differ- 
ences will arise in the MSSM analysis with respect to the SM case. Our conclusion will 
be in fact that the coefficient of the SUSY Sudakov logarithm and the next to subleading 
term will exhibit a drastically different dependence on the parameter of the model, so 
that their possible experimental identification might lead to quite valuable information. 
In particular, we shall reconsider, in this new spirit, the possibility of a determination of 
tan/3, that is the only SUSY parameter on which the coefficient of the SUSY Sudakov 
(linear) logarithm depends, as already shown in previous papers [10,11]. Also, we shall 
show that it is possible to investigate separately the dependence on the next to subleading 
term on the remaining SUSY parameters in the light scenario that we have assumed. In 
fact, we shall also show in some detail what is the range of SUSY parameters that may 
be considered " light" , from the point of view of a Sudakov expansion like ours, in the one 
TeV region. 

We should anticipate at this point the reason why we insist in calling the remaining, 
non logarithmic component of the asymptotic one-loop expansion "next to subleading" 
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term. In fact, we shall show that this remaining component cannot be rigorously con- 
sidered as a constant in the investigated TeV region. To reproduce completely the exact 
calculation, one must add to the constant component an extra non logarithmic energy 
dependent quantity. As we shall show, though, this extra component is, at the expected 
level of accuracy, fairly " small" . The consequence of its presence will be fully taken into 
account and will generate a realistic error in the determination of the SUSY parameter 
that we shall pursue. 

As a first process to be examined in this spirit, we have chosen that of charged Higgs 
pair production. The main reason for this choice is that, from the point of view of the 
involved parameters, this is the simplest process to be considered in the MSSM. Our goal 
is that of moving in future papers to more complicated processes, following a logical chain 
that introduces gradually new parameters not already derived in this effective way. In 
order to provide a self-consistent and rigorous calculation device, we have completed a 
one loop code that contains all the relevant diagrams, in the approximation of treating 
all fermions as massless, with the exception of the third generation quarks. We have ver- 
ified that the complete code does reproduce the correct (known) asymptotic logarithmic 
(DL and SL) behaviour. This complete code has been used to derive an effective next-to 
subleading Sudakov expansion, with which it has been imposed to it to agree to the few 
permille level. The code, that has been named SESAMO (Supersymmetric Effective Su- 
dakov Asymptotic MOde), is already available for use [12]. 

Technically speaking, this paper will be organized as follows: Sec. II will contain 
the relevant asymptotic one loop expansions of the process. Sec. Ill will contain the 
comparison of the complete one loop calculation with the proposed effective fit as a 
function of the SUSY parameters. In Sec. IV, the determination of tan f3 and the study 
of the effect of the remaining parameters on the next-to subleading term will be examined 
and exhibited. A final discussion of our results in Sec. V will conclude the paper. 

II. CALCULATION OF THE PROCESS AT ONE LOOP 

A complete description of the scattering amplitude of the considered process at one 
loop requires the calculation of several classes of diagrams. To make our treatment as 
self-consistent as possible, we shall follow the notations used in a previous reference [10] 
and write: 

A{e + e- -> H + H~) = A Boin (e + e~ -> H + H~) + A c - l \e+ e~ -> H + H~) 

+A s - e -{e + e- -> H + H~) + A in (e + e- -> H + H~) 
+A fin (e + e- -> H + H~) + A hox {e + e~ -> H + H~) + A QED (2.1) 

It is convenient to normalize the general amplitude in the following way: 
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A = % v(e + )({>)(a L P L + a R P R )u(e-) (2.2) 



q 2 



where P LjR = (1 =F 7 5 )/2 and is the outgoing H momentum, so that one writes the 
Born terms as: 

°A - 1 - Aq 2 r 2 V9e\ (2.3) 

with g eL = 2s 2 w - 1, y ei ? = 2s 2 w and 77 = g2 / M . . 

The remaining quantities represent the one-loop perturbative modifications of the 
tree level expression. More precisely, A c - L is the contribution from the usual "counter- 
terms" that cancel all the ultraviolet divergences of the process. In our chosen on-shell 
renormalization scheme, in which the inputs are o:q_ed(0), M w and M z , they are given by 
proper gauge boson self-energies, computed at the corresponding physical masses. Their 
explicit expressions are known (see e.g. [13]), and to save space we shall note write them 
here. A s - e - contains the various internal self-energy corrections; A m ' f m describe the initial 
and final vertex modifications and A box the box contributions. Their "fine" structure is 
summarized in Appendix, where the various components of the separate terms are listed. 
The related complete set of Feynman diagrams is too large (more than 200 diagrams) 
to be drawn here; it can be found e.g. in a recent paper [14]. In our treatment, we 
have discarded those components that give contributions that vanish with the initial 
lepton mass, which reduces somehow our numerical calculations; (actually, we treated all 
fermions as massless, with the exception of the third family quarks). In our approach, 
ultraviolet divergences are produced both by self-energies and by vertices that include all 
external self-energies (we followed the definition proposed by Sirlin [15]). We checked that 
all the ultraviolet divergences in the scattering amplitude are mutually canceling. Finally, 
A® ED represents our choice of the electromagnetic component. In this preliminary paper, 
we were mostly interested in the "genuine" electroweak SUSY contributions at very high 
energies. For this reason, we treated all those virtual contributions with photons that 
would generate infrared divergences by introducing an "effective" fictitious photon mass 
M 7 = Mz- With this choice, our A® ED will contain in fact the difference between these 
"effective" terms and the conventional (massless photon) ones (with the usual request 
of adding the effects of real photon radiation when computing the cross section). This 
gives a contribution to the observables that, for our specific purposes, can be considered 
as "known" since it does not involve any SUSY component, and will therefore not be 
included in our code at least at the moment. Its addition would not represent a problem, 
but would not add anything to the specific investigation of this paper, that is rather 
devoted to the "unknown" component of the scattering amplitude. 

In this spirit, we shall now write the overall logarithmic contributions to the scattering 
amplitude that arise at one loop in a proper configuration of "asymptotic" energy. These 
are by definition the leading terms of an expansion made in a region where all the relevant 
(external and internal) masses are sufficiently smaller than the cm. energy in a way that 
we shall try to make more quantitative in the following Section. In full generality, such 
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logarithmic terms can be of two different origins. The first ones are linear logarithms 
of renormalization group (RG) origin, generated by gauge boson self-energies and repre- 
senting the "running" of the gauge coupling constants. Their expressions are known, and 
we shall write them explicitely but separately from the remaining terms. They do not 
contain any SUSY parameter, but must be carefully taken into account in our approach. 
The second ones are the genuine electroweak logarithms, nowadays generally called "of 
Sudakov type". They can be of quadratic and of linear kind; the quadratic ones come 
from vertices with single W, Z, 7 exchange and from boxes with two W exchange and 
do not contain any SUSY parameter; the linear ones come from the remaining vertices 
and boxes and contain SUSY contributions only from vertices. A very special feature of 
the super symmetric model that we have investigated and of the process that we have con- 
sidered is that the only SUSY parameter that "effectively" appears in the various linear 
logarithms is tan f3, of Yukawa origin, produced by the final vertices with (6, t) exchange, 
that depends on the specific combination [m 2 cot 2 (3 + m\ tan 2 0\ (SUSY mass parameters 
Mi, that could enter other vertex diagrams, would appear in the form log(g 2 /M 2 ), 
is the cm. energy, but with a suitable change of scale they can always be shifted into a 
constant term, as we shall show and discuss). 

In Appendix we list the various logarithmic contributions coming from different dia- 
grams. The convention that we have followed is that of keeping Mw as the scale of the 
double log 2 q 2 coming from two W boxes and single W vertices. For all the remaining 
logarithms we have chosen a common scale Mz, with the exception of the Yukawa vertex 
where the most natural scale appears to us to be that of the top mass m t . This choice 
is arbitrary, it is only dictated by our personal taste and feelings. The consequence of 
this choice will be that of fixing the numerical value of the "next-to subleading" term. A 
pragmatic attitude would be that of verifying that, with this choice, both the logarithmic 
and the next-to subleading terms remain acceptably "small" at the one loop level. This, 
as we shall show, will happen in fact and with this "a posteriori" justification we shall 
retain our choice of scales. 

From a glance to the Appendix one also sees that the linear log q 2 belong to two 
separate classes: the so-called "universal" and the "not universal" ones. The latter are 
produced by boxes and depend on the cm. scattering angle 1?. We have already discussed 
in a previous paper [10] their different relevance at the one loop level, and we shall 
not insist here on these features. We can, though, recall the observed fact [10] that 
these d dependent non universal terms are the only ones that contribute at logarithmic 
level the forward-backward asymmetries at one loop. Since they do not contain any 
SUSY parameter, as we said, there will be no relevant logarithmic contribution containing 
SUSY parameters to the forward-backward asymmetry of H + H~ production. The same 
conclusion can be derived for the longitudinal polarization asymmetry that is known to 
provide only information on the initial state as shown in the previous reference [10]. For 
this reason, we shall concentrate our numerical analysis on the total cross section of the 
process. 

After this rather technical presentation, we are now ready to present the concrete 
numerical analysis. With this aim, we shall now summarize our previous discussion giving 



6 



the asymptotic relative effect on the cross section of the process, by writing it in the form: 

^-Born+lloop ^.Born 

A(« 2 ) = " gB „„ (2-4) 

where in - Born + llo °p W e are retaining only the genuine one loop terms 0(a/n) and not the 
second order terms coming from the square of the one loop contributions since these mix 
with the genuine two loop contributions. 

The logarithmic expansion of A has been derived analytically and is given by the 
expression 

A(, 2 ) = -(^)(^f)log 2 ^ 



2ns 2 w JK l+4s 4 w J ta 

- ( 4^^ )(m ' C0t ^ + m ' tan ^ )l0g ^ 

+ [ 3ttsW w )[ 1+4^ jl ° g Ml 

+A rcm (g 2 ) (2.5) 

where the fourth line contains all single logarithms with the exception of those of Yukawa 
origin (third line). The last term A rem (g 2 ) is the difference between the full one loop result 
and its asymptotic Sudakov expansion including all the double and single logarithms, and 
will be called the next-to subleading term. 



III. VALIDITY OF THE SUDAKOV EXPANSION 

The aim of this Section is that of investigating whether there exists a region of energy 
and of parameters where the rigorous calculation at one loop can be reproduced by the 
effective Sudakov expansion Eq.(2.5), and to determine the relevant features of the next-to 
subleading term A rem (g 2 ). In order to avoid confusion, we anticipate that our analysis will 
be divided into two different sectors. In the first one, summarized in this Section, we shall 
study the dependence of A rcm (g 2 ) on the cm. energy for given values of the parameters 
of the chosen MSSM model. In the second one, we shall study the dependence of A rem (g 2 ) 
on the MSSM parameters for a fixed energy chosen at the representative 1 TeV value. To 
proceed with our analysis, we must now define the MSSM parameters that we shall use 
as input of the calculations. We retained the following five free parameters: 

tan/5, /I, M A , M 2 , M s (3.1) 

i.e. the ratio of the two vevs, the Higgs bilinear coupling, the CP odd Higgs mass, the 
universal gaugino mass and the universal sfermion mass. For this preliminary analysis, we 
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allowed ourselves the simplifications of using the GUT relation Mi = | tan 2 $]yM 2 and of 
setting the trilinear couplings A u = = 0. We computed the Higgs spectrum with the 
code FeynHiggsFast [16] and obtained the masses of charginos, neutralinos and sfermions 
by numerical diagonalization of their mixing matrices. We retained sfermion mixing only 
in the case of the third generation. 

For the purposes of this paper, we have chosen to work in an energy region between 800 
GeV and 1 TeV, considering 1 TeV as an ambitious conceivable final goal of the future LC. 
In our study, we shall assume that a number of precise measurements can be performed 
in that energy range, and we will show what could be the theoretical implications for the 
model in a particular region of its couplings and of its masses. Clearly, all the obtained 
informations and bounds on the various mass parameters could be easily rescaled if the 
analysis were performed in a lower (e.g. 600-800 GeV) energy interval. 

The first problem that we address is that of determining the range of massive MSSM 
parameters for which the Sudakov expansion for A(g 2 ) reproduces the rigorous calculation 
with a simple and understandable expression for the next-to sub leading term A rem (g 2 ). 
For the latter, the simplest possibility would be provided by a constant, and from our 
previous SM experience [10], we would be prepared to the appearance in this case of a 
relatively "large" (i.e. compared to the logarithmic component) quantity. We cannot 
exclude, though, a priori the necessity of including an extra, energy dependent term, that 
should vanish at very large \^q I , but could be numerically relevant and complicated in 
the considered energy range. If this turned out to be the case, the practical validity of 
the Sudakov expansion would be unavoidably reduced, since every tentative effective fit 
to the data would be complicated and affected by large errors, and we shall return on this 
statement in the description of the derivation of tan/3, to be shown in the forthcoming 
Section. 

Given the fact that we have to deal with four massive parameters, we have performed 
four different analyses, in each one of which three parameters were fixed at values that 
we considered "light" with respect to the chosen energy range, and one parameter was 
allowed to vary. In all these analyses we fixed tan/3 at the value tan/3 = 20, that can be 
considered as an average value in the range that we have explored, roughly 2 < tan (3 ^ 40 
(we shall discuss later on the reasons of this choice of the upper value). Fixing different 
values of tan/3 does not change the results of the four analyses, that we chose in the 
following way: 

a) variable /i, fixed M A = 200 GeV, M 2 = 100 GeV, M s = 350 GeV; 

b) variable M A , fixed [i = 300 GeV, M 2 = 100 GeV, M s = 350 GeV; 

c) variable M 2 , fixed fi = 300 GeV, M A = 200 GeV, M s = 350 GeV; 

d) variable M s , fixed /i = 300 GeV, M A = 200 GeV, M 2 = 100 GeV; 

The choices of the fixed values of the parameters in the four analyses are dictated 
by practical reasons and could be reasonably varied. The relatively small value of M 2 
corresponds to the request of having relatively light charginos and neutralinos, although 
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higher values 300 GeV) for M 2 would still be acceptable, as shown by Fig. (2). For 
Ma, the same considerations apply in order to avoid resonance formation. Ms can vary 
in a larger range without apparent problems, and therefore it was fixed at a relatively 
higher value. 

We are now ready to discuss the results of our analyses. This was performed by 
computing numerically the quantity A(g 2 ) with the specific code (SESAMO) that we 
have built for our purposes and which is, we repeat, available for use. From the computed 
quantity we then subtracted all the logarithms of Eq.(2.5) and obtained the remaining, 
next-to subleading term A rem (g 2 ). 

An important preliminary remark must be made at this point. The final goal of our 
analysis is that of showing that a complete one-loop expression can be reproduced by 
a (relatively simple) Sudakov expansion. But a necessary condition to make this search 
useful is also that the one-loop expansion is reliable. To guarantee this condition, we shall 
have to verify that the complete expression of A(g 2 ) (i.e. the full one-loop correction) 
remains acceptably small in the investigated energy region, with e.g. an upper bound 
that we could fix quantitatively at the relative ten percent value. Only after this check, 
our next analysis could be considered as meaningful. 

We are now ready to show the results of the four considered cases, which are the 
following ones: 

1) Case (a): variable fi 

We allowed fi to vary from an initial value of 300 GeV to final values of approximately 
400 GeV (below 300 GeV, we encountered problems in the determination of the sfermion 
mass eigenstates whose discussion seems to us beyond the purposes of this preliminary 
analysis). The full one-loop effect, computed at the representative value ^fq I = 1 TeV, 
remains always very small (below two percent). The values of A rem (g 2 ) in the interval 
800 GeV < v^g 2 < 1 TeV are shown in Fig. (3). One sees that the term is changing by an 
amount roughly equal to five percent of its central value when moving from the beginning 
to the end of the interval. Numerically, this corresponds to a two permille effect that can be 
considered as "fairly small" under our expected experimental conditions of a one percent 
accuracy. When fi becomes larger than roughly 380 GeV, the simplicity of A rem (g 2 ) is 
lost and a complicated energy dependence appears which is due to a resonance effect: 
when we increase fi, one of the two charginos and two neutralinos become progressively 
heavier with masses ~ fi. When fi ~ 400 GeV, we begin to see a kinematical threshold at 
~ 800 GeV that produces a bump in A rem well visible in the figure. Of course, the bump 
shifts to the right as fi is further increases. One observes also from Fig. (3) that A rem (g 2 ) 
is scarcely affected by variations in fi (about half percent when fi varies from 300 to 400 
GeV) and, as a consequence, it does not appear to be a promising candidate for testing 
virtual MSSM effects generated by this specific parameter. 

We repeated our analysis changing only the sign of fi. The obtained curves are prac- 
tically (in agreement with our previous statement) identical with those corresponding to 
positive fi values. For this reason, we shall restrict ourselves from now on to considering 
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conventionally the fj, > scenario. 

2) Case (b): variable Ma 

We have varied M A starting from an initial value of 200 GeV up to approximately 
400 GeV. The full (negative) one-loop effect at 1 TeV remains systematically below a 
relative two percent. As Ma increases, the basic change in the spectrum of SUSY parti- 
cles is that both H + and H° become heavier with approximately M A o ~ M H o ~ M H +. 
This means that in the plots we have to take into account the kinematical constraint 
v^2 > 2M H+ . In Fig. (4) we show the behaviour of A rem for 200 GeV <M A < 400 GeV. 
Once again, we notice that in the considered energy interval A rcm (g 2 ) remains "essen- 
tially" constant, with relative extreme variations of five percent (or less) from its central 
value. This simple pattern would be lost for larger Ma values > 400 GeV, due to the 
aforementioned Higgs production threshold. One notices in this case that the dependence 
of A rcm (g 2 ) on Ma is sizable: to a variation of Ma from 200 to 300 GeV there corresponds 
a variation in A rcm (g 2 ) of almost two percent, that would be visible at the proposed LC. 

3) Case (c): variable M 2 

Here, we vary M 2 from 100 GeV to 400 GeV. The full one-loop effect at 1 TeV is again 
very small (below two percent) and negative. In the considered range of variation of M 2 , 
there are one chargino and one neutralino with masses increasing approximately as M 2 
and reaching the 400 GeV value at which some resonance structure can be observed in 
the plots of A rcm shown in Fig. (5). It can be noticed that, for M 2 larger than about 360 
GeV, the shape of this function is definitely not constant with energy and clear bumps 
can be seen, due to resonance effects associated to the heavy gauginos. As in the case of 
/x, we observe that A rem (g 2 ) has very small sensitivity to the variations of M 2 (about two 
permille variation for a 100 GeV shift in M 2 ). 

4) Case (d): variable M s 

We have varied Ms from 100 GeV to about 400 GeV. The full one-loop effect is again 
very small and negative (below ~ 3 — 4 percent). The values of A rcm (g 2 ) shown in Fig. (6) 
are again essentially (up to a relative 5 percent) constant in the considered energy inter- 
val. One sees that A rem (g 2 ) is sensitive to M s : 100 GeV of variation in this parameter 
corresponds to ~ 1.5 percent variation in A rem (g 2 ). 

To summarize, we have considered the range of validity of a logarithmic Sudakov 
expansion with an extra next-to subleading term which is "essentially" constant in a 
certain domain of the massive MSSM parameters, for cm. energies in the range 800 GeV— 
1 TeV. In our approach, to be "essentially" constant means to be well approximated by 
the central value, with a few permille error. We have verified that this requirement is met 
for values of all the parameters below, approximately, 350 GeV. This gives a quantitative 
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illustration of how "light" the MSSM parameters should be in order that a typically 
asymptotic expansion like the Sudakov one that we investigate might hold at energies in 
the 1 TeV range. We can remark that, as matter of fact, the validity of the expansion 
corresponds to values of the masses meeting the naive request Mf/q 2 < 1CT 1 . This 
remains true even if several resonances appear close to y/q z /2. In other words the energy 
dependence of this process appears to be reasonably flat in the considered domain. We also 
conclude that the values of the "essentially" constant components are relatively "large" 
and of opposite sign (positive) with respect to the (negative) logarithmic contribution 
(roughly, they are of comparable size, although always smaller than the logarithmic ones). 
As a result of this cancellation the overall one-loop effect is sensibly reduced, becoming 
systematically of few percent size. This situation reproduces exactly the one that we met 
in the analogous study of the SM case, even if the number of free parameters in the MSSM 
is so much larger. 

A final question that we asked ourselves was that of whether the approximate con- 
stancy of A rem (g 2 ) could be due to the relative smallness of the investigated energy domain. 
To answer this question, we have extended our analysis to very large ^fq I values up to 
10 TeV, for an illustrative fixed set of MSSM parameters 

tan/3 = 30, fj, = 300, M A = 200, M 2 = 100, M s = 350, (3.2) 

The results are shown in Fig. (7). One sees that for energies beyond ~ 2 TeV the 
values of A rem (g 2 ) remain practically (i.e. to less than one permille) constant. Although 
we already know that at such energies the one-loop approximation is probably not valid, 
we consider this result as a check of the asymptotic validity of the expansion (and also of 
the numerical code that we have used). 

In conclusion, we have seen that in the 1 TeV range a 1-loop Sudakov expansion 
with an "essentially" constant term is valid for the defined "light" SUSY scenario. The 
fact that A rcm (g 2 ) is not rigorously constant is expected to produce therefore reasonably 
"small" effects, that we shall try to evidentiate in full detail in the forthcoming Section 4. 

IV. STUDY OF THE MSSM PARAMETERS 
A. Determination of tan/3 

The first question that we address is the relevance of the Sudakov expansion for the 
problem of determining tan f3 under the assumption that the mass scales are in the "safe" 
range that we have just discussed. In other words, A rcm is only approximately constant 
and we want to understand quantitatively how much this can affect an attempt to de- 
termine tan /3 from the logarithmic slope of the cross-section, in which it appears in the 
combination shown in Eq.(2.5). 

With this aim, we begin by subtracting explicitly from A all the "known" logarithms, 
i.e. all the terms in the Sudakov expansion with the exception of the Yukawa contribution. 
Therefore, we define the quantity 
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A(g 2 ) = F(tan/5) log q 2 + A rcm (g 2 ), (4.1) 
F(tan/5) = - a (m 2 cot 2 (3 + m 2 tan 2 /5). 

If, in a definite scenario, the shape of A rem (g 2 ) turned out to be flat, then it would be 
conceivable to approximate it by a term that is constant with respect to q 2 . In such a 
simple case, we could try to fit the measured values of the residual effect A(g 2 ) with a 
logarithmic expansion in q 2 of the form 

A fit logg 2 + £ fit . (4.2) 

The result of the fit, A& t , can be compared with F at the value of tan/5 we are working. 
The difference 5F = A &t — F is an error in the estimate of F that has two components: 
SF = S stSut F + 5 sys F. The first term, 5 stat F is simply due to the fact that we assume 
a certain finite experimental precision on each measurement. The second term, <5 sys is 
the most important and is a systematic error due to the fact that A rem is not constant 
with respect to the energy. For instance, if A rem were exactly energy independent, we 
would find S sys F = 0. The error 5F can be converted into an error on the estimate of 
the interesting parameter tan/5. If 5F/F is enough narrow to allow a linearized analysis, 
then we have simply 



SF _ tan/5 F' <5tan/5 
~F ~ F tan/5 

or 

5tan/5 _ 1 tan 4 /5 + (m t /m b ) 2 5F 
tan/5 ~~ 2 tan 4 /5- (m t /m b ) 2 ~F 



(4.3) 



(4.4) 



The zero in the denominator corresponds to the value tan/3 = \Jm t / m& ~ 6.2 at which 
the function F attains its maximum and the sensitivity to tan /5 is the smallest due to the 
flatness of F. Notice also that for tan/5 beyond 15 — 20, we have 

5F 1 5tan/5 . . 

"f ~ 2 15T/T (45) 

In the following discussion we shall analyze in a quantitative way the feasibility of such 
a procedure in the framework of specific scenarios. With this aim, we have assumed the 
existence of 10 equally spaced experimental measurements in the range 800 GeV — 1 TeV 
with a relative 1 percent precision and have generated them by means of our numerical 
code (if only iV points are available, all the numerical results concerning the statistical 
component of the error on tan/5 must be increased by a factor yJl0/N). 

L: Very Light SUSY 

(L) variable tan/5 /j, = 300, M A = 250, M 2 = 100, M s = 350. (4.6) 
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The full effect at f TeV is given in Fig. (8). showing that it remains below a fO percent 
for tan/? ^ 40. The curves for A rem are given in Fig. (9) showing that A rcm (g 2 ) depends 
effectively on tan/3, remaining "essentially" constant in the considered energy range. The 
plot of the relative error in the determination of tan (3 is shown in Fig. (10). As we said, 
the extra error bars are due to the fact that A is determined with a statistical error due 
to the assumed 1% accuracy in the cross section measurements. One sees that the main 
source of error is actually due to the departure of A rem (g 2 ) from its constant value. 
We have subsequently considered two more scenarios, defined as: 

A: Light SUSY 

Here, we increase the masses in the gaugino sector. 

(A) variable tan/5 fi = 300, M A = 250, M 2 = 200, M s = 350. (4.7) 
B: Light SUSY with larger \x 

(B) variable tan/5 fj, = 400, M A = 250, M 2 = 200, M s = 350. (4.8) 

In Fig. (11) we combine the results for the various scenarios. In the figure, we have 
also shown the vertical lines corresponding to the safe perturbative bound corresponding 
to a 10% full one loop effect. One can see that, for tan/3 larger than 20, a determination 
of this parameter to better than a relative 40 percent would be possible. For values larger 
than ~ 30, the error would be reduced below a remarkable ~ 10 percent limit. This 
would remain true even in the worst case of relatively heavy SUSY scenario (3) with 
e.g. Ma = 250 GeV, and would represent to our knowledge a valuable possibility of 
determining this fundamental MSSM parameter in the region of high values where it is 
known [17] that accurate measurements are rather difficult. 

B. Visible effect of the remaining parameters 

Our logical scheme for extracting information from charged Higgs production would 
now proceed in the following way. Once the proposed determination of tan f3 from mea- 
surements of the slope of the cross section were completed, we would return to the remain- 
ing term A rcm (g 2 ) and estimate the effects on it of the remaining parameters assuming 
a precise measurement at a fixed energy, typically 1 TeV. A preliminary request will be 
that of taking into account the error on the previous determination of tan (3. Following 
the illustrations of Section 3, we shall optimistically assume that tan/5 has been deter- 
mined at a "convenient" value, i.e. one where the relative error is of the ten percent 
size. For purposes of illustration, we shall chose tan (3 = 30 ± 3 from now on. Using this 
value as a given input, we can examine which information on the remaining parameters 
can be obtained from the determination of A rem . This determination will be affected by 
two sources: a purely experimental one from the measurement at 1 TeV, treated under 
the usual assumptions, and the input error on tan/5 measured from the slope. It is not 
difficult to see (e.g. looking at Fig. (9) and considering Eq.(2.5)) that the latter will affect 
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the determination of A rem by a tolerably small (few permille) error. We shall take it into 
account in what follows within qualitative limits, not to make this indicative treatment 
too involved. 

The plan of our forthcoming study has been remarkably helped by the observation that 
we already made, which shows that, in practice, A rem remains "essentially" unaffected by 
variations of fi and of M2 in the considered " light" scenario. This simplifies our approach, 
reducing it to the "essential" parameters, that are Ma and M5. We have thus fixed [i, M2 
at conventional values (/i = 400 GeV, M 2 = 100 GeV) and drawn the contour and the 
surface plots in the (M s , M A ) variables shown by Figs. 12,13. 

A few, necessarily qualitative, comments are now appropriate e.g. from a glance to 
Fig. 12. The various curves correspond to variations of A rcm at 1 TeV. The spacing between 
two curves is a shift in A rem of 5 permille, which corresponds roughly to one half or our 
expected error on this quantity and, in our Figure, defines a certain bi-dimensional "tube" 
whose slope and width depend on the parameter domain and fix the corresponding domain 
bounds on (M A , M s ). One notices that, independently of the value of A rem , there would 
be a kind of orthogonal situation. For small 4 percent) A rem values, one would feel 
the effect of Ms with a certain accuracy (of about 50 GeV) without practical effect from 
Ma- For larger A rem values, the opposite situation would appear, and effects of Ma could 
be felt to the previous (about 50 GeV) accuracy. These accuracies are certainly much 
worse than the expected precisions on Ma, Ms from direct production (roughly, a relative 
1-2 percent). However, in our opinion, these curves could still be rather meaningful for 
a possible non trivial consistency test of the model. Assuming in fact that both Ma and 
Ms have been determined in a range between 200 and 350 GeV with a precision of, say, 
five GeV, the point in the (Ma, Ms) plane that corresponds to these values must lie on 
the "correct" curve that corresponds to the measured value of A rcm . In Fig. 12 we have 
drawn for illustration purposes two points that correspond to typical couples of "light" 
values M A = 250 GeV, M s = 330 GeV and M A = 330 GeV, M s = 250 GeV with the 
corresponding assumed experimental error. One sees that a measurement of A to the 
relative one percent accuracy would be of scarce use for Ma = 330 GeV, M s = 250 GeV, 
but would provide a quite stringent test of the model for the symmetrical couple of 
values. Thus, depending on the experimental results on these masses, the relevance and 
the motivations of the previous analysis at one TeV might become definitely enhanced. 

V. CONCLUSIONS 

The main conclusions that may be drawn from our analysis of the charged Higgs pro- 
duction process are, in our opinion, the following 

1) For this process, in the ~ 1 TeV energy region, an effective one- loop description 
with a Sudakov expansion implemented by an "essentially" constant next-to sub leading 
term reproduces the rigorous calculation in a "light" SUSY scenario where all the relevant 
mass parameters of the process are roughly below the common 350 GeV value. The over- 
all one-loop effect remains systematically under control (below a safe few percent limit) 
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in this region and seems to provide a reliable description of the process. 

2) A satisfactory determination of tan f3 from an accurate measurement of the slope of 
the cross section in a region close to and below 1 TeV would be possible for large values 
(> 20) of tan/3, with an error that takes realistically into account the "small" deviations 
of the next-to subleading term from a constant value. 

3) The next-to subleading term "essentially" depends, once given tan/3 from the mea- 
sured slope, only on the two mass parameters Ma, M s . Depending on the measured values 
of the parameters, this could provide a simple but rather stringent test of the MSSM. 

4) For the purposes of a test of the MSSM, the charged Higgs production process ex- 
hibits special simplicity features that make it, in our opinion, a very promising candidate. 
This property would suggest to extend our study to cases of non minimal SUSY models, 
in particular models whose Higgs couplings and structure are either different or richer. 
We have completed our examination using the dedicated code SESAMO; our next step 
will be that of generalizing our analysis to the similar case of neutral Higgs production. 
This will require an extension of our code, on which work is already in progress. 
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APPENDIX A: LIST OF CONTRIBUTIONS AND ASYMPTOTIC 

EXPRESSIONS. 

In this Appendix we give the list of the various one loop diagrams that have to be 
retained for a practical computation (for example we discard all diagrams that contribute 
proportionally to the light lepton and quark masses). We follow the decomposition given 
in eq.(2.1). 

Gauge boson self-energies 

These are the standard and supersymmetric bubbles and seagull diagrams involv- 
ing gauge bosons (7, Z, W), goldstones, ghosts, Higgses, fermions, charginos, neu- 
tralinos and sfermions. They contribute the quantities A cd '(e + e~ — > H + H~) and 
A s - e ie + e~ -> H+H~), as explained in refs. [13,14]. 

Initial vertices 

The diagrams contributing A m (e + e~ — > H + H~) are vertices with three internal lines 
sketched in Fig. (1) and external self-energies. The list of vertices (a, b, c) is: (e, 7, e), 

(e, Z, e), (u, W, v), (W, v, W), (e, X °, e), fc, X + , *e), (x°, e, *e, X + )- 

final vertices 

The diagrams contributing A* m (e + e~ — > H + H~) are vertices sketched in Fig. (1) and 
external if ± bubbles as well as seagull diagrams involving the gauge boson-gauge boson- 
scalar-scalar couplings. The list of vertices (a,b,c) is: (H, 7, H), (H, Z, H), (H, W, H), 
(W, H, W), (/, /', /), ( X , X , X), (H, H, Z), (H, H, H), (/, /', /), where // and X 
represent either charged or neutral states. The list of seagull diagrams (a, b) is (H + ^) , 
(H+Z), (H°W), (h°W), (A°W). 

boxes 



The contributions to A box (e + e — > H + H ) arc box diagrams denoted clockwise by 
starting from the line running between e~ and e + according to Fig. (1). The list of 
boxes (a,b,c,d) is: (vWH°W), (e7^+7), (eZH + Z), (eZH+^f), (e-fH + Z), (ux + X°X + ), 
(ex°X + X°), (x°e~Lve~ L ), (x~i>e L i>), {yx + X°e L ), (e L X°X + v)- 

Asymptotic expressions 

The complete expressions have been included in the code SESAMO. Below we only 
give the results involving leading (DL) and subleading (SL) logarithms. 

Using the normalizations defined in Eq.(2.2), we can write: 
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7T 

where <5a^ c -, <5a™ and £af n are the one-loop corrections to a^ orn . 

The asymptotic contributions from the intermediate 7, Z self-energies are: 

SaT - l ~]fM 2sAw log g 2 54 e - - log g 2 

those from initial e + e~ lines: 



2 log, 2 -log 2 '' 



2 



and those from final H + H~ lines and boxes: 

- T^T f2 log q 2 - log 2 ^1 + f2 log g 2 - log 2 " 



64s^ w V b " b Ml) 32s* w c* w V b " b 
3 

(m 2 cot 2 /3 + m\ tan 2 /3) log g 2 



1 j 1 — cos-# j 2 1 j 1 — cos-i? j 2 
16sfyCw ° S 1 + cos-tf ° S 2 



- 4 ( 2 ,og « 2 - log2 A) + isU ( 2 log?2 - log2 A, 

m 2 cot 2 /3 + m 2 tan 2 /3) logg 2 



1 1 — cos-i? 2 
log • log q 



Actv 1 + cosd 
where $ is the cm. angle between initial e~ and final H~ momenta. 
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FIG. 4. Variable M A : plot of A rcm for various M A . The bottom line has M A = 200 GeV. 
The other curves have increasing M A by steps of 20 GeV. 
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FIG. 7. Behaviour of A rem in the ultra high energy regime. 
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FIG. 8. Variable tan (3 in the (L) scenario: full effect at 1 TeV 
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FIG. 9. Variable tan/3 in the (L) scenario: A rem . The bottom line has tan/3 = 10, the other 
curves have increasing tan/3 by steps Atan/3 = 2. 
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FIG. 12. Contour plot of 100A rem (ireV) in the plane (M A ,M S ). The contour lines corre- 
spond to values between and 8.5 increasing from left to right by steps of 0.5. Also shown 
are the boxes corresponding to the points (Ma, Ms) = (250 ± 5 GeV, 330 ± 5 GeV) and 
(M A , M s ) = (330 ± 5 GeV, 250 ± 5 GeV). 
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FIG. 13. Surface plot of 100A rcm at 1 TeV in the plane (M A ,M S ). 
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